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ABSTRACT 

We quantify the level of foreground contamination in the QMAP Cosmic Microwave Background 
(CMB) data with two objectives: (a) measuring the level to which the QMAP power spectrum mea- 
surements need to be corrected for foregrounds and (b) using this data set to further refine current 
foreground models. We cross-correlate the QMAP data with a variety of foreground templates. The 30 
GHz Ka-band data is found to be significantly correlated with the Haslam 408 MHz and Reich and Reich 
1420 MHz synchrotron maps, but not with the Diffuse Infrared Background Experiment (DIRBE) 240, 
140 and 100/zm maps or the Wisconsin H-Alpha Mapper (WHAM) survey. The 40 GHz Q-band has no 
significant template correlations. We discuss the constraints that this places on synchrotron, free-free 
and dust emission. We also reanalyze the foreground-cleaned Ka-band data and find that the two band 
power measurements are lowered by 2.3% and 1.3%, respectively. 

Subject headings: cosmic microwave background - diffuse radiation - radiation mechanisms: thermal 
and non-thermal - methods: data analysis 



1. INTRODUCTION 

Quantifying Galactic emission in a cosmic microwave 
background (CMB) map is interesting for two different 
reasons. On one hand, the CMB is known to be a gold 
mine of information about cosmological parameters. Tak- 
ing full advantage of this requires accurate modeling and 
subtraction of Galactic foreground contamination. On the 
other hand, the high fidelity maps being produced as part 
of the current CMB gold rush offer a unique opportunity 
for secondary non-CMB science. This includes a greatly 
improved understanding of Galactic emission processes be- 
tween 10 and 10 3 GHz. 

This paper is motivated by both of these reasons. The 
QMAP experiment (Devlin et al. 1998; Herbig et al. 1998; 
de Oliveira-Costa et al. 1998b, hereafter dOC98b) is one 
of the CMB experiments that has produced a sky map 
with accurately modeled noise properties, lending itself to 
a cross-correlation analysis with a variety of foreground 
templates. We present such an analysis in §^ and §|[ 
then compute the corresponding correction to the pub- 
lished QMAP power spectrum measurements in §^ and fin- 
ish by discussing the implications for Galactic foreground 
modeling in ^ 

2. METHOD 

The multi-component fitting method that we use was 
presented in detail in de Oliveira-Costa et al. 1999 (here- 
after dOC99), so we review it only briefly here. The joint 
QMAP map from both flights consists of N = 3164 (Ka- 
band, 26 to 36 GHz) and 4875 (Q-band, 36-46 GHz) mea- 
sured sky temperatures (pixels) j/j. We model this map 
as a sum of CMB fluctuations Xi, detector noise rij and 
M Galactic components whose spatial distributions are 
traced in part by external foreground templates. Writing 



these contributions as A-dimensional vectors, we obtain 



y = xa 



n, 



(1) 



where X is an N x M matrix whose rows contain the 
various foreground templates convolved with the QMAP 
beam (i.e., would be the i th observation if the sky had 
looked like the j th foreground template), and a is a vector 
of size M that gives the levels at which these foreground 
templates are present in the QMAP data. 

We treat n and x as uncorrelated random vectors with 
zero mean and the X matrix as constant, so the data co- 
variance matrix is given by 



(yy T >-(y)(y T > = (xx r ) + (nn r ), 



where 



E 

1=2 



2£+l 

47T 



(2) 



(3) 



is the CMB covariance matrix and (nn T ) is the QMAP 
noise covariance matrix (a detailed description of the cal- 
culation of (nn T ) is presented in dOC98b). We use a 
flat power spectrum Ci cx + 1) normalized to a 

Q flat = (5CV47T) 1 / 2 = 30 ^K (dOC98b). We model 
the QMAP beam as a Fisher function with FWHM= 
V8 ln2cr = 0.9° for Ka-band and 0.6° for Q-band, which 
gives a We « e- - 2 ^ 1 )/ 2 . 

Since our goal is to measure a, both x and n act as 
unwanted noise in equation ([!]). Minimizing \ 2 = (y — 
Xa) T C _1 (y — Xa) yields the minimum- variance estimate 
of a, 

a = [X T C" 1 X] _1 X T C 
with covariance matrix 

£ = (a 2 ) - (a) 2 = [X T C- 1 



X 



(4) 



(5) 
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Fig. 1. — Top: from left to right, the 100/xm, Has, H a and 
QMAP Ka-band maps are shown, each in their own units. They 
are centered roughly on the North Celestial Pole. Bottom: from 
left to right, the three templates above are multiplied by the a- 
values from Table 1 and summed. These four maps and the CMB 
map all use the fi¥L scale shown by the color bar, so the the best 
fit map of the foreground contribution (lower right) is seen to be 
quite subdominant to the CMB signal (upper right). To make the 
CMB plot less noisy, the Wiener filter described in dOC98b has been 
applied to all eight maps before plotting. 

The error bars on individual correlations are therefore 

1/2 

Aft, = . This includes the effect of chance alignments 
between the CMB and the various template maps, since 
the CMB anisotropy term is incorporated in (xx T ). 

3. BASIC RESULTS 

We cross-correlate the QMAP data with two different 
synchrotron templates: the 408 MHz survey (Haslam et 
al. 1982) and the 1420 MHz survey (Reich 1982; Reich 
and Reich 1986), hereafter Has and R&R, respectively. 
To study dust and/or free- free emission, we cross-correlate 
the QMAP data with three Diffuse Infrared Background 
Experiment (DIRBE) sky maps at wavelengths 100, 140 
and 240/im (Boggess et al. 1992) and with the Wisconsin 
H- Alpha Mapper (WHAM)Q survey (Haffner et al. 1999). 
For definiteness, we use the DIRBE 100/mi channel when 
placing limits below since it is the least noisy of the three 
DIRBE channels. Three of our templates are shown to- 
gether with the Ka-band QMAP map in Figure 0. Most 
of our interesting results come from the QMAP Ka-band, 
since the Q-band was substantially noisier (the opposite 
was true for the Saskatoon experiment - see de Oliveira- 
Costa et al. 1997, hereafter dOC97). 

Before calculating the correlations, we convolve the tem- 
plate maps with the QMAP beam function. We also re- 
move the monopole and dipole from both the templates 
and the QMAP maps. As a consequence, our results de- 
pend predominantly on the small scale intensity variations 
in the templates and are insensitive to the zero levels of 
the QMAP data and the templates. 

Table 1 shows the coefficients a and the corresponding 
fluctuations in antenna temperature in the QMAP data 
(AT = aacai, where acai is the standard deviation of 
the template map). Statistically significant (> 2a) cor- 
relations are listed in boldface. Note that the fits are 
done jointly for M — 3 templates. The DIRBE, Haslam 
and H a correlations listed in Table 1 correspond to joint 



100/im— Has— H a fits, whereas the R&R numbers corre- 
spond to a joint 100/im— R&R— H a fit. Only the two syn- 
chrotron templates are found to be correlated with the 
Ka-band, while no correlations are found for the Q-band. 
Repeating the analysis done for two different Galactic cuts 
(20° and 30°) indicates that the bulk of this contamination 
is at latitudes lower than 30°. 

As in dOC97, de Oliveira-Costa et al. 1998a (hereafter 
dOC98a) and dOC99, the cross-correlation software was 
tested by analyzing constrained realizations of CMB and 
QMAP instrument noise. From 1000 realizations, we re- 
covered unbiased estimates a with a variance in excellent 
agreement with equation (ph. As an additional test, we 
computed x 2 = (y — Xa)"C _1 (y — Xa) and obtained 
X 2 /N w 1 in all cases. Including a synchrotron template 
lowered x 2 by a significant amount (18 for R&R and 9 
for Has), whereas adding the other templates resulted in 
insignificant reductions A% 2 ~ 1. 

4. IMPLICATIONS FOR CMB 

The lower right map in Figure ^ shows Xa T , i.e., our 
best fit estimate of the foreground contribution to the 
QMAP Ka-band. To quantify the foreground contribution 
to the published QMAP power spectrum measurements, 
we repeat the exact same analysis described in dOC98b 
after subtracting out this map, i.e., with the map y re- 
placed by y — Xa T . 

The dOC98b band powers were computed by expanding 
the map in signal-to-noise (S/N) eigenmodes (Bond 1995; 
Bunn and Sugiyama 1995), weight vectors that solve 

Table 1 

Correlations for Ka- and Q-band data. 



b & v 


Maps (a) 


a ± AaW 




4li 

Gat 


AT [/tK]M 


P 


All Sky 


100/im 


-3.5 ± 6.9 


-0.5 


3.1 


-10.8±21.4 




Ka-band 


Has 


6.6 ± 2.6 


2.5 


4.3 


28.4±11.2 


-2.8 




R&R 


0.3 ± 0.1 


2.6 


141.4 


42.4±14.1 


-2.6 




H a 


14.6 ± 9.5 


1.5 


1.0 


14.6± 9.5 




All Sky 


100/im 


-1.7 ± 5.4 


-0.3 


3.9 


-6.7±21.4 




Q-band 


Has 


0.7 ± 2.8 


0.2 


6.1 


4.2±17.2 


-2.7 




R&R 


-0.05± 0.09 


-0.6 


182.3 


-9.4±16.3 


-3.0 




H a 


-6.0 ±11.2 


-0.5 


1.1 


-6.3±11.8 




b > 20° 


Has 


6.1 ± 2.7 


2.3 


2.2 


13.9± 6.2 


-2.8 


Ka-band 


R&R 


0.3 ± 0.1 


2.2 


82.5 


24.0±10.9 


-2.7 


6 > 30° 


Has 


4.1 ± 5.2 


0.8 


1.4 


5.8± 7.4 


-2.9 


Ka-band 


R&R 


0.4 ± 0.3 


1.2 


28.6 


10.5± 8.9 


-2.6 



(a) The DIRBE and Haslam correlations listed in this table 
correspond to joint 100/tm— Has— H a fits, whereas the R&R 
numbers correspond to a joint 100/im— R&R— H a fit; 

(b) a has units /iK(MJy/sr) _1 for the 100/im template, /tK/K 
for the Has template, /tK/mK for the R&R template and /iK/R 
for the H a template; 

' c ' cFGai has units of the templates; 
< d > AT= (a± Aa)a Ga i. 



1 Details also at 



http://www. astro, wise, edu/wham/ . 
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the generalized eigenvalue equation 

(xx T )b, ; = A i (nn T )b, : 



(6) 



When the b's are sorted by decreasing eigenvalue A, they 
tend to probe from larger to smaller angular scales. The 
S/N expansion coefficients are shown in Figure |^, and are 
seen to be substantially smaller for the foregrounds than 
for the CMB. As described in dOC98b, we obtain a sta- 
tistically independent power estimate from the square of 
each mode and then average these individual estimates 
with inverse-variance weighting to obtain the band power 
estimates listed in Table 2 and shown in Figure || The Ka- 
band (30 GHz) band powers are seen to drop by less than 
a few percent when the foreground signal is subtracted, 
whereas the Q-band (40 GHz) contamination is too small 
to quantify. The Ka-band correction is slightly smaller on 
small angular scales: 1.3% instead of 2.3%. This is ex- 
pected, since diffuse Galactic foregrounds are expected to 
have a redder power spectrum than the CMB. 

As a side benefit, our statistically significant detection 
of foregrounds allowed an independent confirmation of the 
QMAP pointing solution. We reanalyzed the QMAP data 
set with the pointing solution offset by ±0.5° in azimuth, 
and found the highest correlations for the original point- 
ing. 

5. IMPLICATIONS FOR FOREGROUND MODELING 

Figure || shows a compilation of measured correlations 
between foreground templates and CMB maps at various 
frequencies. Comparisons are complicated by the fact that 
the foreground level AT = aacai depends strongly on 
galactic latitude via acai- We therefore plot a instead, 
i.e., the factor giving the frequency dependence of emission 
per unit foreground. Such measurements were used to nor- 
malize recent foreground models such as those of Bouchet 
& Gispert (1999) and Tegmark et al. (1999). Below we 
discuss how our QMAP results affect such models. 

5.1. Synchrotron 

Writing the frequency dependence as a oc v" and re- 
calling that the correlation coefficients are, by definition, 
a = 1K//xK= 10 6 for Has at 408 MHz and a = lmK//iK= 
10 3 for R&R at 1420 MHz, we obtain the spectral index 
limits —2.7 < < —3.3 for the Ka-Has correlation and 
-2.6 < (3 < -2.8 for the Ka-R&R correlation. These val- 
ues are slightly steeper than the canonical sub-GHz slope 
of -2.7 < (3 < -2.9 (Davies et al. 1998; Platania et 
al. 1998), but consistent with a steepening of the spec- 
trum of cosmic ray electrons at higher energies (Rybicki 
and Lightman 1979). 



Table 2 

The angular power spectrum before and after foreground removal. 

The band powers ST e = [£(£+ l)C e /2n] 1 / 2 have uncorrelated 
errors and window functions whose mean and rms width are given 
by £ ef f and A£ just as in dOC98b. 
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ST after 
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127 


55 


59+ 6 7 


58±? 




40 60 
Mode number 
Fig. 2. — S/N eigenmode coefficients of the Ka-band data (top) 
and combined foregrounds (bottom). The wiggly curves show the 
a running average of the rms fluctuations with bin length 5 and 
25 The smooth curves show the theoretically expected rms for flat 
power spectra with Qji a t =0, 10, 20, 30 and 40 fi¥L. 
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Fig. 3. — Power spectrum from QMAP experiment. The fore- 
ground removed points (in magenta) is shown with the original 
points (in black). 



The relatively high QMAP synchrotron signal seen in 
Figure |^ could be interpreted as slight spatial variability 
of the frequency dependence (Tegmark 1999), but may also 
have other explanations. For instance, the worst striping 
problems in the Haslam map are right around the North 
Celestial Pole, which may have caused Saskatoon to un- 
derestimate the true synchrotron level there (dOC97). 

5.2. Spinning dust & free- free emission 

An important question is whether the DIRBE-correlated 
signal seen by so many experiments (see Figure ^, top) 
is due to dust-correlated free-free emission (Kogut et 
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Fig. 4. — Galactic contaminants as function of frequency. Fre- 
quency dependence of DIRBE-correlated emission (top), Haslam- 
correlated emission (middle), and H a -correlated emission (bottom). 
From left to right, the dust points are from Tenerife 10 GHz and 
15 GHz (dOC99: filled triangles; Mukherjee et. al. 2000 Table 2 
b > 20°: open triangles), 19 GHz (dOC98a), this work (the upper 
limit), DMR 31 GHz (Kogut et al. 1996), Saskatoon (dOC97), DMR 
53 and 90 GHz (Kogut et. al. 1996). The Synchrotron points are 
from the same experiments as the dust points at the corresponding 
frequencies. The entry in the bottom plot are from this work. Up- 
per limits are 2-<r. The dashed model curves show the frequency 
dependence characterizing (spinning + vibrating) dust (top), syn- 
chrotron radiation (middle) and free- free emission (bottom). The 
normalization of the synchrotron model a oc v~ & is fixed by the fact 
that, by definition, a = lK//iK= 10 6 at 408 MHz, and the dust 
model must similarly reproduce DIRBE at 100/tim. 



eq. (|5|) shows that we can interpret S _1 as the covari- 
ance between the various templates with dimensionless 
correlation coefficients = S~ 1 (S^ 1 ST. ) -0 ' 5 . Like in 
dOC99, the DIRBE maps were found to be almost per- 
fectly correlated, and essentially uncorrelated (r 2 < 3%) 
with the radio maps. The Has and R&R maps are cor- 
related with r ~ 83% for b > 20°. As a new result, we 
obtain a marginal correlation r « 0.2 between the DIRBE 
maps and the H a template. Since the statistical prop- 
erties of these maps are not accurately known, we com- 
puted error bars by repeating the analysis with one of 
the templates replaced by 2 x 2 x 72 = 288 transformed 
maps, rotated around the Galactic axis by multiples of 
5° and/or flipped vertically and/or horizontally. The ac- 
tual correlation was found to be larger than 85% of these, 
showing that the correlation is not significant at the 2a 
level: a = (0.25 ± 0.19) R/MJysr" 1 (la). This result is 
significantly smaller than that recently found by Lagachc 
et al. (2000) for their DIRBE-WHAM correlation done in 
a different region of the sky, but compatible with other 
marginal dust-_ff Q correlations (McCullough 1997; Kogut 
1997). 

This poor correlation is a challenge for the pure free- 
free hypothesis, which maintains that microwave emission 
traces dust because dust traces free-free emission. A cross- 
correlation analysis with large frequency and sky coverage 
will hopefully be able to unambiguously determine the rel- 
ative levels of free-free and dust emission in the near fu- 
ture. 

We would like to thank Matias Zaldarriaga for encourag- 
ing us to complete it. Support for this work was provided 
by the Packard Foundation, NSF grants PHY 92-22952 
& PHY 96-00015, and NASA grant NAG5-6034. We ac- 
knowledge the NASA office of Space Sciences, the COBE 
flight team, and all those who helped process and analyze 
the DIRBE data. 



al. 1996) or spinning dust grains (Draine & Lazarian 1998). 
The turndown at low frequencies suggests a spinning dust 
interpretation (dOC99), but an analysis using improved 
Tenerife data (Mukherjee et al. 2000) may have re-opened 
this question^. Our present results cannot settle the issue, 
but offer a few additional clues. 

If free-free emission is responsible for the correlation, 
then substantial H a emission would be expected as well. 
Figure |] (bottom) shows the expected correlation a for the 
case of 8, 000K gas (Bennett et al. 1992). It is seen that the 
iJ Q -corrclation, although marginal at best, is consistent 
with the theoretical curve. Moreover, Table 1 shows that 
the limit on the dust-correlated emission is low enough to 
be compatible with a free- free origin (i.e., with the ~ 15/iK 
iJa-correlated component). 

To clarify this issue, we computed the correlation be- 
tween the dust and H a maps. As described in dOC99, 

2 As seen in Figure M, the Mukherjee et al. results for b > 20° 
also shown a turndown Detween 15 and 10 GHz. Figure 7 in their 
paper looks different because it is an average including data down 
to b = 0. This may still be consistent with spinning dust being 
the worst foreground for the (Galaxy cut) MAP data, since free-free 
emission is believed to be more concentrated in the Galactic plane 
than dust. 
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